BAR4R-CONF-04/50 
SLAC-PUB-10646 
hep-ex/0000000 
July 2004 



A search for the decay K^uu 

The BABM Collaboration 
February 7, 2008 



Abstract 

We present a search for the rare flavor-changing neutral-current decay B'^ K^vV based on 
data sample of 82fb~^ collected with the BABAR detector at the PEP-II B-factory. Signal events 
are selected by examining the properties of the system recoiling against either a reconstructed 
hadronic or semileptonic charged B decay. Using these two independent samples, a combined 
limit of B{B~^ — > K'^uV) < 5.2 x 10~^ is obtained at the 90% confidence level. In addition, by 
modifying the particle identification criteria, a limit of B{B^ tt^uu) < 1.0 x 10""* obtained using 
the hadronic B reconstruction. All results are preliminary. 
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Flavor-changing neutral-current transitions such as 5 ^ svV and b — > dvV occur in the Standard 
Model (SM) via one-loop box or electroweak penguin diagrams with virtual heavy particles in 
the loops. Therefore they are expected to be highly suppressed. Because heavy non-SM particles 
could contribute additional loop diagrams, various new physics scenarios can potentially lead to 
significant enhancements in the observed rates [T]. Theoretical uncertainties on 6 ^ svV are much 
smaller than the corresponding b — > si~^i~ modes due to the absence of a photonic penguin contri- 
bution and hadronic long distance effects The SM — > K^uV branching fraction has been 
estimated to be (S.SlgJ) ^ 10"^ [21 E], while the most stringent published experimental limit is 
B{B^ K^uV) < 2.4 X 10~^ at the 90% confidence level (C.L.) PJ. b ^ duV processes are ad- 
ditionally suppressed relative to 6 — >• svV by the ratio of the Cabibbo-Kobayashi-Maskawa matrix 

elements \Vtd?/\yts? U- 

In this work we report the results of a search for the exclusive decay mode — s- K~^uV. By 
modifying the particle identification (PID) criteria used in the search, we additionally obtain a 
limit on the related decay B^ tt^vV. 

The data used in this analysis were collected with the BABMi detector ^ at the PEP-II 
asymmetric-energy e'^e~ storage ring. The results are based on a data sample of 88.9 million 
BB events, corresponding to an integrated luminosity of 82 fb~^ collected at the T(4£') resonance. 
An additional sample of 9.6 fb~^ was collected at a center-of-mass (CM) energy approximately 
40 MeV below BB threshold which is used to study continuum events, e^e~ qq {q = u, d, s and 
c). Charged particle tracking and dE/dx measurement for particle identification (PID) are provided 
by a five-layer double-sided silicon vertex tracker and a 40-layer drift chamber contained within 
the magnetic field of a 1.5— T superconducting solenoid. Charged K—tt PID separation of greater 
than 3(7, over the momentum range of interest for this analysis, is provided by a ring-imaging 
Cherenkov detector (DIRC). The energies of neutral particles are measured by an electromagnetic 
calorimeter (EMC) consisting of 6580 CsI(Tl) crystals. The magnetic flux return of the solenoid is 
instrumented with resistive plate chambers in order to provide muon identification. A full BABAR 
detector Monte Carlo (MC) simulation based on GEANT4 [7] is used to evaluate signal efficiencies 
and to identify and study background sources. Charge conjugate modes are implied throughout 
this paper and all kinematic quantities are expressed in the CM frame (i.e. the T(4S') rest frame) 
unless otherwise specified. 

The presence of two neutrinos in the final state precludes the direct reconstruction of the 
B^ K^vV signal mode. Instead, the B~ meson from an T{AS) B^B~ event is reconstructed 
in one of many semileptonic or hadronic decay modes, then all remaining charged and neutral 
particles in that event are examined under the assumption that they are attributable to the decay 
of the accompanying B. 

The B~ reconstruction proceeds by combining a candidate with either a single identified 
charged lepton or a combination, X^^^, of charged and neutral hadrons. The resulting semileptonic 
and hadronic charged B samples are referred to as B"^ and B^^^ throughout this paper, 
candidates are reconstructed by selecting combinations of identified pions and kaons which yield 
an invariant mass within approximately 3cr of the expected mass in the modes K~7r~^, K~tt~^tt^ 
and K~7r~^7r~^Tr~ . For -B^ad reconstruction, KgTr~^7r~ is also used. 

Photon candidates are obtained from EMC clusters with laboratory-frame energy greater than 
30 MeV and no associated charged track. Photon pairs which combine to yield 77 invariant mass 
between 115 MeV/c^ and 150 MeV/c^ and total energy greater than 200 MeV are considered to be 
tt" candidates. 

B~^ candidates are reconstructed by combining a candidate having a momentum p^o > 
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0.5 GeV/c with a lepton candidate of momentum > 1.35 GeV/c that satisfies either electron or 
muon identification criteria. The invariant mass, m£)£, of the D^i candidate is required to be greater 
than 3.0 GeV/c^. candidates are selected using the quantity cos 9 b^di, which represents the 
cosine of the angle between the inferred direction of the reconstructed B~y and that of the lepton- 
combination, described by the four vector, {EuijPDi). Under the assumption that the neutrino 
is the only missing particle, 



2K 



cos tiB,De = 



beam 



• E 



m 



B 



m 



Dl 



2 Ipd^I 



-^beam ^B 



(1) 



where is the nominal B meson mass and -Bbeam and y E^^^^ — m\ are the expected B meson 
energy and momentum, respectively. Combinatorial backgrounds can produce values, \cos6B,Dt\ > 
1. In order to maintain efficiency for B~ — > D*^i~hi decays in which a tt^ or photon has not been 
reconstructed as part of the D(. combination, we retain events in the interval —2.5 <cos 9b, d£< l-l- 
However events are vetoed if a charged vr consistent with a D*+ transition is identified. If more 
than one Di candidate is reconstructed in a given event, the candidate with the smallest \cos 9b, ml 
is retained. 

Reconstructed B^^^ decays are obtained by combining a reconstructed candidate with a 
hadronic system X^^^ composed of up to five mesons {7r^,K^ and tt''), including up to two vr^ 

candidates. We define the kinematic variables tjies = \J -^beam ~ Pb ^^"^ ^ — -Ebeam where 
Pb and Eb are the momentum and the energy of the B^^^ candidate. The system is selected 
by requiring that the resulting B^^^ candidate lies within —1.8 < /S.E < 0.6 GeV. If multiple B^^ 
candidates are identified in an event, only the one with AE closest to zero is retained. The rriES 
distribution of reconstructed -B^ad candidates is shown in Fig. ^p. B^^^ candidates in the signal 
region, 5.272< mss <5.288 GeV/c^, are used for the B^ K'^uV signal selection. Candidates in 
the sideband region, 5.225< ttt-es <5.265 GeV/c^, are retained for background studies. 




D° mass (GeV/c^) 



nrigs (GeV/c^: 



Figure 1: (a) The mass distribution for — s- K^tt^ decays used for B~^ reconstruction. Data 
are shown as points and the total background MC is shown as a solid histogram, (b) The tties 
distribution of B^^^ events for data (points) and BB MC (solid histogram). Continuum background 
has been subtracted from the on-resonance data using off-resonance data and the hatched histogram 
represents the estimated combinatorial background from BB decays. 



Combinatorial backgrounds from continuum events are reduced in both the B~^ and B^^ sam- 
ples by requiring | cos^tI < 0.8, where 9t is the angle between the thrust axes defined by the B~y 
or i?bad daughter particles, and by all other tracks and clusters in the event. Continuum events 
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peak at | cos 9t\ = 1, while the distribution is approximately flat for BB events. Backgrounds from 
QED processes are strongly suppressed by the B~ reconstruction procedures and are negligible in 
this analysis. 

The B~ reconstruction efficiency for B~^ K^vV signal events is determined from signal MC 
simulation after validating the yield from B^B~ MC simulation against data. This procedure 
compensates for differences in the -B^ad reconstruction efficiency in the low-multiplicity environment 
of B^ K^uV events compared with the generic B^B^ environment. The B~~ reconstruction 
efficiency in MC is additionally validated by comparing the yield of events in which a B^ — > D^l^v 
has been reconstructed in addition to the -B^ad -^sT- 

The B~^ and B^^^ reconstruction procedures result in raw yields of approximately 5800 B^^/ fb^^ 
and 2200 B^^^/ fb~^. Relative systematic uncertainties of 4.5% (7%) are estimated for the overall 
^sT (^had) yields. 

Events that contain a reconstructed B are examined for evidence of a B'^ K'^uv decay. 
Tracks and EMC clusters not already utilized for the B~ reconstruction are assumed to be the 
daughters of the signal candidate B decay. Signal candidate events are required to possess exactly 
one additional charged track with charge opposite that of the reconstructed B~ . The track is 
additionally required to have momentum pK > 1.25 GeV/c and to satisfy K PID criteria. 

In addition to this track, i?+ K^uV events contain an average of ~ 200 MeV of EMC energy 
from hadronic shower fragments, photons from unreconstructed D* — > D^j/ir^ transitions in the 
B'^ candidate, and beam-related background photons. The total calorimeter energy attributed to 
the signal decay, -Bextra, is computed by summing all EMC clusters that are not associated either 
with the decay daughters of the B~ or with the signal track. Signal events are required to have 
-E'extra < 250 MeV. The -Eextra distributions are shown in Fig. 121 for B^^ and i^^ad events with one 
additional track which has been identified as a kaon. The B^^ analysis additionally requires that 
there are six or fewer clusters contributing to -Ecxtra, and that no pair of these clusters can be 
combined to form a vr*^ candidate. 

The total B^ K^vV signal selection efficiencies, including the B~ reconstruction, are esti- 
mated to be EK = (0.115 ± 0.009)% for B'^ and er = (0.055 ± 0.005)% for B^^^ events. The quoted 
errors are the quadratic sum of statistical and systematic uncertainties. Theoretical uncertainties 
in the energy spectrum result in a 1.3% uncertainty on the signal efficiency. This uncertainty is 
evaluated by comparing the pK spectrum of B^ — > K^vV MC events generated with a phase-space 
model with the models given in [HI |2j. Additional systematic uncertainties associated with the 
B^ K^uV signal candidate efficiencies include the single track efficiency (1.3%), PID (2%) and 
EMC energy modeling (3.8% for B~y and 2.3% for -Bj^ad)- EMC energy modeling systematic 
is determined by evaluating the effect of varying the MC -Eextra distribution within a range repre- 
senting the observed level of agreement with data in events with a reconstructed B^ D^i'^u (for 
the B~^ sample) and in samples containing two or three additional tracks (for the -B^ad sample). 

Background events can arise either from B^B^ or continuum events in which the B~ candidate 
is constructed from a random combination of particles, or peaking background events in which the 
accompanying (or in the case of B^, at least the D^) has been correctly reconstructed. 

In the B~^ analysis, purely combinatorial backgrounds are estimated by examining sideband 
regions of the reconstructed invariant mass distribution, rrf^°, defined by 2>a < |m^o° — "i_do| < 
10(7 as shown in Fig. ^ for the K~tt'^ mode. The sideband yields are scaled to the signal 

region under the assumption that the combinatorial component is flat throughout the mass 
distribution. This assumption has been validated using samples of events in which two or three 
charged tracks not associated with the B~ reconstruction are present. The total combinatorial 
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background in the analysis is estimated to be > 3.4zbl.2. Although the peaking background 
prediction in the analysis have been studied in MC and are shown in Figs Inland 01 the peaking 
background in the final selection is not subtracted. 
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Figure 2: The -Eextra distribution for K^vV (left) -B^ad (right) B^^ events. Events are 

required to have a reconstructed B~ and exactly one additional track which has been identified as 
a kaon. No other signal selection cuts have been applied. The data and background MC samples 
are represented by the points with error bars and solid histograms, respectively. The dotted line 
indicates the B~^ K~^vV MC prediction with arbitrary normalization. 



In the -B^ad analysis, the combinatorial background can be reliably estimated by extrapolating 
the observed yields in the ttt-es sideband region into the rn-Es signal region, indicated in Fig. 
yielding 2.0 it 0.7 events. The quoted uncertainty is dominated by the sideband data statistics, but 
includes also the uncertainty in the combinatorial background shape which is estimated by varying 
the shape over a range of possible models. The peaking background in the -Bj^ad analysis consists 
only of B~^B~ events in which the -B^ad been correctly reconstructed, and is estimated directly 
from B^B^ MC simulation. MC yields are validated by direct comparison with data in samples of 
events in which the full signal selection is applied, except that either -Eextra > 0.5 GeV, or more than 
one charged track remains after the B~ reconstruction. Uncertainties in the peaking background 
are dominated by the MC statistical uncertainty (42%). Other systematic errors include the overall 
B~ yield (7%), the remaining charged track multiplicity (5%), the particle mis-identification rates 
for the selection (6.3%), and the EMC energy modeling (8%). The total peaking background 
in the -B^ad analysis is estimated to be 1.9 ± 0.9. The total (combinatorial+peaking) background 
in the B^^ analysis is estimated to be A''^^ = 3.9 ± 1.1 events. 

The B~^ K~^vV branching fraction is calculated from 

ATobs _ ji^bg 

B{B+ ^ K+vV) = -^^ (2) 

Nb± ■ ek 

where N'^^ is the total number of observed events in the signal region. Nq± = (88.9 it 1.0) x 10^ 
is the estimated number of B^ mesons in the data sample and ek is the total efficiency. A total of 
N'^^ = 6 (3) B~^ K^uV candidate events are found in data in the B~^ i^had) analysis. The pK 
distributions for B^ K^uV signal events in the B~y and B^^ analysis are shown in Fig. [21 

Branching fraction upper limits are computed using a modified frequentist approach, based 
on which models systematic uncertainties using Gaussian distributions. For both the B~^ and 
■^had searches, B~^ — > K~^vV limits are set at the branching fraction value at which it is estimated 
that 90% of experiments would produce a yield which is greater than the number of signal events 
observed. Limits of B{B+ K+uV)^x < 7.0 x 10"^ and B{B^ K+uV)^^^ < 6.7 x 10"^ are 
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Figure 3: The pK distribution for (a) -B^ad ^^'^ (^) ^sl events after applying the fuh — > K'^uV 
selection except for the pK > 1-25 GeV/c requirement. The data and background MC samples are 
represented by the points and solid histograms respectively. The dotted line indicates the signal 
MC prediction for — > K^uV with arbitrary normalization. 



obtained for the B^^ ™d -^had searches respectively. Since the two tag B samples are statistically 
independent, we can combine the results of the two analyses to derive a limit of B{B^ — > K^i'V) < 
5.2 X 10"^ at the 90% C.L. 

We also report a limit on exclusive B^ tt^uV branching fraction using only the B^^^ sample. 
The same methodology as for the B^ — > K'^vV search is applied to the B^ — > tt'^vV search except 
that the single additional track is required not to satisfy either kaon or electron PID criteria. The 
-E'extra and distributions for B^ tt^uV are shown in Fig.^ The overall B^ 'k^vV selection 
efficiency is estimated to be = (0.065 ± 0.006)%, where the quoted uncertainties include an 
estimated 2% PID uncertainty, and other contributions to the systematic uncertainty are similar to 
B^ K^uv. The peaking and non-peaking backgrounds are estimated to be 15.1 it 3.1 events and 
9.0±1.8 events respectively, with similar systematic uncertainties to the B~^ K^uV analysis. The 
search selects N!^^ = 21 candidates in data with an estimated total background of N^^ = 24.1 ±3.6, 
resulting an upper limit of B{B^ 7r+i/I7)had < 1-0 x 10"'' at the 90% C.L.. 
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Figure 4: The -Eextra (a) and p^^ (b) distributions for B~^ — > tt'^uV in the B^^^ sample. Events shown 
in the -Eextra distribution are required to have a reconstructed B~ and exactly one additional track 
satisfying the pion selection requirements. The distribution has all signal selection requirements 
applied other than the cut. The data and background MC samples are represented by the 
points and the solid histogram respectively. The dotted line indicates the signal MC prediction 
with arbitrary normalization. 



We see no evidence for a signal in either of the reported decay modes. The B{B^ K^vv) 
limit reported here is approximately one order of magnitude above the SM prediction and represents 
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the most stringent experimental limit reported to date. 
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